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coordinates in the parallel space. The numerical data obtained with the two 
methods are in impressive agreement. 
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1 Introduction 



The purpose of this paper is to calculate typical two-loop IR-divergent vertex 
diagrams by two independent methods. One of them is based on general for- 
mulae for asymptotic expansions of Feynman integrals in momenta and masses 

^ (see for a brief review) and subsequent use of conformal mapping and 
summation by Pade approximations [||, |2|. The general simple formulae have 
been proven Q at least for the Feynman integrals off the mass shell? *However 
it is quite natural to expect that the same off-shell formulae hold as well for 
the pure large mass limit when there are no large momenta, as it has been con- 
firmed in ref. This conjecture is likely to hold even for Feynman integrals 
which possess IR and collinear divergences from the very beginning. To check 
the validity of this conjecture we shall calculate IR-divergent diagrams 7 and 8 
shown in Fig. 1 (the labelling is according to Fig. 3) by use of the large 
mass expansion and compare results in a wide range of external momenta with 
results based on numerical integration. 

The project is motivated by the demand to provide the necessary integrals 
for the process Z — > 65. The before-mentioned divergences appear in the limit 
rrih — > 0, which is a natural limit to take since rrih is small compared to the 
other scales in this process. The case of a finite m^, mass can in principle also 
be handled by both methods, see e.g. Ref. for the momentum expansion 
method applied to a finite diagram (Case 5, ms = mg = M, all others zero, 
see Fig. 1). In this case, however, one obtains less precise results because 
of the low thresholds, which are eliminated in the present approach due to the 
factorization of logarithms. For the numerical method the case 7^ presents 
no difficulty. 

The numerical methods, as introduced in |3|, ^, IC, are in the present 



paper extended to the case of IR and collinear divergences. We will explain this 
extension below. The two methods used in this paper are of complementary 
nature. While for the method based on the large mass expansion the equal 
mass case is simpler than the case of different masses, the situation is reversed 
if we apply the two-fold integral representation. The presence of poles of the 
type \\m.mx^mi ^/i'fni — 1^2) in the integral representations results in numerical 
instabilities which have to be handled with due care. Corresponding poles in the 
first method can be avoided from the very beginning by calculating the Taylor 
coefficients analytically in terms of rational numbers and C{2) for the equal 
mass case. In the case of nonequal masses the Taylor coefficients contain such 
poles of high order (see below), which cancel in the limm^^mj- Numerically a 



close approach is possible due to the use of a multiple precision package [12|. 

The idea of the numerical approach is now well documented in the literature. 
We will thus restrict ourselves to describe the necessary changes demanded 
by the presence of IR and collinear divergences. These changes are simple: 
one merely has to find subtraction terms to absorb the divergences, which are 



*Note that for some typical on-shell limits explicit formulae for asymptotic expansions have been 
recently presented and applied in Ref. M. 
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sufficiently easy to allow for a direct evaluation. So there was no new conceptual 
challenge present in the numerical method, but solely a technical hurdle to be 
overcome. 

The paper is organized as follows. First, we discuss the method based on 
the large mass expansion for Cases 7 and 8, considering the equal as well as 
different masses in Case 7. We then discuss the relevant subtraction terms 
for the numerical method. Once they are understood, all the previous cases 
can be easily obtained. We summarize the results in a set of tables, and com- 
pare accuracies and CPU times for both methods. We finally conclude that 
the agreement of data obtained by the two independent methods confirms the 
usefulness of both. 




Pi -P2 
generic 




Case 7 



Case 8 



Figure 1: On-shell infrared divergent planar diagrams with zero thresholds (solid lines 
massive, dashed lines massless). 



2 Case 7 

2.1 Equal masses 

Let Fy{pi,P2, M; e) be the Feynman integral corresponding to Case 7, with 
ms = rui = ms = 0, mi = m2 = me = M, Pi = P2 = 0, Q = Pi — P2- We imply 
that the factor — m"^ + iO) corresponds to a line, l/(27r)'* is not included. 

For convenience we divide our Feynman integrals by i7r'^/^r(l + e)/i~^^ per 
loop, where d = 4 — 2e is the space-time dimension and /x the scale parameter 
of dimensional regularization (although we shall not usually write down ^ 
explicitly). 

Although the diagram has IR poles, up to let us apply the general 

formula for the large mass expansion |^, ^ (as it was applied for similar IR- 
finite planar diagrams in Ref. [^). There are contributions from two subgraphs 
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in this formula: the graph F = 71 itself and 72 = {126} (we denote subgraphs by 
collections of their lines). The first, "naive", contribution is the formal Taylor 
expansion of the initial diagram in external momenta. It reduces to two-loop 
vacuum graphs with numerators, two massive and one massless line. This term 
possesses IR divergences that were not present from the very beginning. 

The second contribution is nothing but the Taylor expansion of the heavy 
triangle (with pi+k2,P2 + k2 as external momenta, ^2 being the loop momentum 
of the light triangle) inserted into the light triangle. Now, this heavy triangle is 
a scalar function of three variables: (pi + A;2)^, (^2 + ^2)^ and q^. When perform- 
ing its Taylor expansion, the factor (pi + A;2)^ leads to a cancellation of one of 
the lines in the light triangle. This produces a one-loop massless diagram with 
its external momentum on the lightcone. It is zero within dimensional regular- 
ization. Thus we come to the conclusion that only terms without such factors 
survive and one can put /c2 to zero in the expansion of the heavy triangle. So, 
the term under consideration happens to be just a product of two factors: the 
light triangle with zero masses and the Taylor expansion of the heavy triangle 
with pi,P2,Pi — P2 as external momenta. The light triangle is equal to 



£2 r(l -2e) (-g2 _io)i- 



(We shall later on omit — iO in —q^ — iO, for brevity.) 

The Taylor expansion of the heavy triangle is also calculated explicitly: 

- ^ ^ V J!i±i±£)_ro2/2M2)" (2) 

The second term does not involve UV divergences, in contrast to the large 
mass expansions considered in Ref. [^. In large mass expansions (as well as 
large momentum expansions off the mass shell) induced IR divergences are 
usually cancelled by induced UV divergences. Here we have an example of a new 
cancellation: the induced IR divergences are cancelled by collinear divergences 
that are present in the second term. To a large extent this happens when we 
put to zero the above mentioned one-loop massless diagram on the lightcone 
where the UV divergences are cancelled by the collinear divergences. 

Thus the contribution of the second subgraph takes the form (with the 
normalization described above) 

^ ^ ' '-Vm')-^(mV/x2)-^ ^ r(i-£)2 



qiM-i^ ^ ' e'^T{l + e)T{l-2e) 

r(n + l + e) 



00 



;^Q2(n + l)(2n + l)!!^'' / > 

If an initial quantity is finite the cancellation of poles serves as a good check 
of the asymptotic expansion. In our case, it is reasonable to check that the poles 
in the sum of the terms in our expansion are the same as in the initial diagram. 
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To see that the poles in e in the sum of the two contributions are the same 
as in the initial Feynman integral let us apply the following expression for the 
pole part in e of the product of three propagators, namely, l/k"^, l/{k'^ — 2pik) 
and — 2p2k) considered as a distribution in k: 



7r2- 



— { (-^ + 7E^) S^'^Hk) + i £ ^[5{k - zpi) + 5{k - zp2) - 25ik)] 



(4) 

It is easy to observe that the term with 6^'^\k) exactly corresponds to the 
pole part of the second contribution. Thus, the term with the integral over z 
should give the pole part of the naive contribution. Let us therefore calculate 
the integral over z of the rest of the diagram (which is nothing but a one-loop 
triangle diagram with the external momenta (1 — z)pi,p2 — zpi, plus a similar 
contribution, with p2 instead of pi). Representing the result in the form of a 
Taylor expansion in the external momenta we conclude that the pole part in e 
of the naive part is equal to 

where s„ = Si{n) and we introduce 

n I 

Sk{n) = j:-,. (6) 

This property has served as a check in our calculation. 

The double pole in e comes only from 72 . Using (|3[) we obtain the following 
explicit expression for the coefficient at 



2\2 



arcsin"' Wg2/4M2 . (7) 



The naive part is calculated as described in see also |14]. The two-loop 
bubble integrals with equal masses and one zero mass are again of the type 
as used in |2|, i.e. expressible in terms of T-functions as a result of which 
the Taylor coefficients of the naive part are essentially rational numbers (apart 
from a ({2) as for the second contribution. These Taylor coefficients have 



been calculated with FORM [15| and the Fade approximants in turn by means 
of REDUCE 1 16], which allows an easy change to floating point numbers of 
arbitrary precision. Adding up the finite parts, the results are presented in 
Table 1 with the following normalization. In addition to — 7r'^r(l -|- e)^//"^*", we 
extract the factor r(i-2£) choose fi = M. 

2.2 Two different masses 

Froceeding as before in Case 7 with two different masses, 777,3 = itT'A = = 
0, 7771 = ^772 = -^1) TTT'S = ^2, wc come to the following result for the contri- 
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bution of the second subgraph: 



^[Mf) {2n-^ + l)\ (M| - Mf)!^^ / ^'^ ^ " 

Calculation of the pole and finite parts allows for the representation 

^ + ^^^^+ro-nL + ir,L^ (9) 

£^ £ 2 

where L = In M| + ln(— g^) and 



-j^ oo 



-0 = 5Mf?|:"'hVMf)" 



2n+l 

3z*+i(i-n)!(2n -« + !)! 

+3s„+it„+l + ln^(l - 2) + - S^+i - -Un+l " "iSi-nU-n)} , 

-I 00 f n—1 I ■ -I \ 1 

1 V-^U„2/,^2Nn V- (n-2-1)! 



2z*+i(i-n)!(2n-i+ 1)! 



y — TT7 ^^^77 77 ( " ln^(l " 2) " 2 ln(l - Z)Si-n 



-, 00 



X 



(n-i-1)! 2n+i (^^^^ _ ^ _ ^^^^^ I 

z^+i(2n-i + l)! z^+^{i-n)\{2n-i + iy. J' 



where = C(2) - 5*2 (n) and = -2C(3) + Ssin), Skin) given by Eq. (|) 
and 2 = 1 — MI/M^. The importance of presenting the above formulae lies in 
their easy numerical evaluation. If instead we would use the expanded form, 
30 coefficients could not even be compiled anymore simultaneously. 

We also calculate (in the same way as in the case with equal masses) a quan- 
tity, the pole part of which is equal to the pole part of the naive contribution: 
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1 1 ^ 



y ^ „=i j=i J 

f2n-j + i + r{j -i-n-l + e) (M|)"+^-J'+i-^ 
' ^ 2n-j + l j (J^^^)! (M| - M2)2n-i+i+2 

^"^+^ A- + A (-l)T(n - j - i + g) (Mf )-+^-— - I , 
+ ^ 1^ J J (2n-j-i + l)! (M| - Mf p+^+i J ^'^ ^ • ^'"^ 

Explicitly, this gives 

12" 

n=l 



j 1^ t:;^V 2n-j + l / (i-i)!(n- j + z + l)!z2n-i+*+2 

"yVi + A /n + i\ Injl - z) + Sj \ 

t^o V ;(2n-i-z + l)b^+^+i ^\ j )il{n+l-i)lz^+^+^ j ^ 

Furthermore this representation of the pole part has been used as a check 
of the naive part. 

Finally there remains to calculate the finite contribution of the naive part. 
Again we proceed as in |^] and |jl^ . The difficulty arising now is a complication 
in the bubble integrals with two different masses and one zero mass (e.g., = 
0): 



VB[a,p,-f,mi,m2,m3) = 



(mf)2 J {kf-mlnkl-min{h-k2r-mlp-' 

(12) 

Apart from F-functions these contain now the hypergeometric function 



2Fi{a,b,c; z) with the argument z = 1 j < 1 [17| (z = and 2-^1 = 1 

rrij — 

in the equal mass case!). Therefore it does not seem to be advisable anymore 
to use FORM for the evaluation of the Taylor coefficients since too complicated 
expressions would arise. Instead we extended a program developed for arbi- 
trary nonzero masses M based on the multiple precision FORTRAN by D.H. 



Bailey [12|. In this framework the occurrence of a zero mass, i.e. the additional 
IR divergences cause the following complication: while for arbitrary nonzero 
masses only the bubble integrals 

Vb{1, I,n,mi,m2,m3) = F(l, 1, n, mi, 7712, m3)+- — „ ^, n > 2, 

£ [n — l)(n — 2) (mg)!.""^^ 

(13) 
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(F(- • •) being the finite part) are divergent like ^ (apart from Vb{1,1,1) and 
Vb{1, 1, 2) which have also ^ -terms), in the present case all the needed bubble 
integrals VB{l,m,n) with m > 1 and n > 3 have a ^-term, which is to be 
calculated by means of the hypergeometric function 2Fi{a, b, c; z) with integer 
indices a, b,c {d = 4) and c = a + b + 1. This 2F1 can be expressed in terms of 
ln{l — z) and powers of z and 1 — z. 

Knowing the divergent parts of the bubble integrals explicitly, the recursion 
for their finite part can be performed and the complete contribution to the 
finite part of the diagram can be calculated. 

Finally two further remarks are in order. The first concerns the evaluation 
of the second subgraph. As one observes from the Taylor series of ro,ri and 
r2 in (9), these contain high powers of In the equal mass case {z = 0), the 
expansion of the ln{l — z) to high powers cancels all inverse powers of z and the 
finite coefficients as described in Sect. 2.1 are obtained. Numerically, however, 
due to the use of the multiple precision program of Bailey |12| for mass ratios 
close to 1 (e.g. z ~ 10~^) still stable results in close agreement with the equal 
mass case can be obtained. 

The second remark concerns adding up of the various contributions. There 
are two possibilities: one can sum the different series ( for ro,ri etc. ) by 
means of Fade 's, multiply the various results with the corresponding powers 
of L (see e.g. (9)) and sum — or one can add up on the level of the coefficients, 
multiplying the Taylor coefficients with the powers of L, and then apply Fade 's. 
We used the latter procedure since the coefficients are known to high precision 
(of the order of 50 decimals) . The summed series were of course by far not that 
precise and if cancellations occur, one looses more precision than necessary. 
That in this manner correct results are obtained is demonstrated by comparing 
the results with those by the numerical method (see tables). 

While in the present paper adding different Taylor series with kinematical 
factors is performed for the calculation of one diagram only, the same procedure 
will also work when adding scalar amplitudes from various Feynman diagrams 
- and may indeed work in an optimal way. 

For the pole parts we do not show tables since they are either given explicitly 
in terms of known functions or arbitrarily many coefficients can be obtained 
easily and thus arbitrarily high precision. The results for the finite part are 
presented in Table 2. We use the same normalization as in the case with the 
equal masses and choose ^ = Mi . 



3 Case 8 

Let here i*r(pi,p2 5 M; e) be the Feynman integral corresponding to Case 8, 
with mj = 0, i = 1, . . . ,5, mg = M. Now, in its large mass expansion, we 
have contributions from four subgraphs: 71 = r,72 = {3456}, 73 = {126}, 
and 74 = {6}. In the last two contributions, the factor (||) is again naturally 
factorized. Straightforward calculation leads to the following result for the sum 
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of them: 

_ 1 1 i r(i-£)3 

~ (M2)i+£ (-^2)i+s^r(l-2e) 

The calculation of the contribution of 72 is similar to the corresponding 
calculations for Cases 1 and 5 in ref. |p, with the following result: 

1 11°° 

jM^sJZ^e-Y./n\e){e/Mr. (15) 

where 

c(f)(e) = ^"Y^\_i){h+i2+n,)/2 (n-(ii + i2-n3)/2)! 

H,i2,n3>0, n+i2+n3 even j3>0 ~ + «2 + n-i)/2)\ 

*1 + ^2 "'"'"'3 

^l!^2!^(^l +^2 - j3)^(n -^2 + j3)^(-^i +^2 +^3) 

((ng - h)/2)\{{h + 12- i3)/2)!((ii - 22 + j3)/2)!((-ii +12+ i3)/2)! 

r(i-£)r((n + i2-j3)/2 + i-£) 

r(l + e)r((ii + i2 - i3)/2 + 2 - 2e) 
xC7(l + i2, 3 + n + (ii - ^2 + n3)/2; {h + ^2 + n3)/2), (16) 

and 



r(ri + r 2 - 5 - f )r(s - r2 + 

r(ri)r(. + f) 



C^(n,r2; .) = ' '^./.J"^; ' ' ^^ (17) 



Note that this expression is obtained from the corresponding contribution of 
Case 5 in by the change 

C(2 + ii+i2,2 + n- iii+i2- n^)/2; (n + ^2 + "-3)/2) 
^ C(l + ^2, 3 + n + (ii - i2 + n3)/2; (ii +12 + n3)/2). 

To see that the poles in e in the sum of all the four contributions are the 
same as in the initial Feynman integral let us once again apply eq. (^). Now 
we observe that the term with 5^'^'' (k) corresponds to the pole part of the sum 
of the contributions from 73 and 74 given by (|l^). Thus, the term with the 
integral over z should give the pole part of the sum of contributions from 71 
and 72- So, we perform integration over z of the one- loop triangle diagram with 
one non-zero and two zero masses. Furthermore we apply general formulae for 
the large mass expansion to this very triangle and eventually conclude that the 
pole part in £ of the sum of the contributions from 71 and 72 should be equal 
to the pole part of the following quantity: 

g2(M2)i+^.r(l + e);^^"-|^^^^ J J(^/^) 

I o 1 1 r(l-£) ^ SnTjn+l-e) 2/^2^n (.c,^ 

^ I^^^^^h r(n + 2-2.) 1^ ^ ■ ^^^^ 
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Note that double poles in ([l^) cancel so that there must be as well a cancellation 
of double poles in the sum of the contributions of 71 and 72. Again this has 
been checked and also that the single pole part agrees with the one obtained 
by direct calculation of the sum of 71 and 72. 

For Case 8, the double pole in e originates only from the sum of the contri- 
butions of 73 and 74. From (|l^ we have the following formula for its coefficient: 

^ {ln(-gVM=^)ln(l + +Li2(-gVM')} • (19) 

Concerning the finite part, the situation is similar as in Sect. 2.1, Case 7, 
with equal masses. The bubble integrals are of the same type as in Case 5 
and the calculation of the Taylor coefficients have been done with FORM and 
the Fade 's with REDUCE. Again we present no data for the divergent parts 
since arbitrary precision can easily be obtained from the expansion. The results 
for the finite part are presented in Table 3 where the same normalization as in 
Case 7 is used. 



4 The Numerical Method 



The two- fold integral representation derived in |18, y, 1C| cannot be naively 
applied to cases infected by genuine IR or collinear divergences. Nevertheless, 
all these divergences can be handled by appropriate subtraction terms. One 
observes a few new characteristics in such cases: 

• the divergences are most easily absorbed using lightcone coordinates for 
internal momenta in the parallel space, 

• some of the remaining domains of integration are unbounded, cf. Fig. (5). 
We will comment on these features in the next three subsections. 



4.1 P reparat ions 

The graph we have to calculate is 



2d-12 



where C{k) denotes the inner one-loop triangle graph with external momenta 
shifted by k (Fig. |). 

For our calculation, we choose the rest frame of the decaying particle and the 
outgoing particles moving along the x axis. This is the natural reference frame 
for parallel-/orthogonal space splitting. Together with the condition that both 
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outgoing momenta are on the hghtcone, the external momenta can be expressed 
through a single parameter e (cf. Fig. ||): 

q = qi+q2 = (2e,0,0) 
qi = (e, e, 0) 

92 = (e,-e,0) (21) 
As usual we parameterize the loop momenta / and k as follows: 

l'' = {k,h,l±) and kf" = {ko,ki,k^), (22) 

further we define 

8 = 11 and t = k'i. (23) 

In contrast to |l^ we apply a more symmetric substitution to linearize the 
propagators in Iq, li, ko and ki (the Jacobian gives a factor 2 for each loop 
momentum) : 

^0 lo + h 
h h — Iq 
ko ko + ki 

ki ^ ki-ko, (24) 

which is equivalent to expressing the propagators in lightcone variables from 
the beginning. After this substitution the propagators become 

Pi = Alo{li + e) - s - mj + iO 

P2 = All ik-e) - s-ml + iO 

P3 = a(^Iq -\- kQ){li + ki) - 8 -t - ^/sViz - ml + iO 

Pa = ikoiki - e) - t - ml + iO 

P5 = 4:ki (ko + e) - t-ml + iO 

Pq = Akoki -t-ml + iO (25) 

Since we are interested in the collinear divergent case, we now set m4 and 772.5 
equal to 0, mg should be kept small for the finite contribution as a regulator. 
The limit will be made later analytically. The calculation is valid for arbitrary 
masses mi, 7712 and 7773, except for the trivial case 7771 = 777-2 = '^s = 0. 

In this representation the collinear divergence along the lines A;o = and 
/ci = for t = is obvious. To cure it, we adopt the following subtraction 
scheme for C{k) = C{kQ, ki,t): 

C{koM.t) = C(A;o,A;i,t)-C(fco,0,0)-C(0,A:i,0) + C(0,0,0) (/) 
+C(fco,0,0)-C(0,0,0) (//) 
+C(0,A;i,0) -C(0,0,0) (I/I) 
+C(0,0,0). {IV) 

(26) 
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Contribution (I) will be finite and can be calculated in d = 4 dimensions, 
since we subtracted out the coUinear divergences and added again the twice 
subtracted infrared divergence at k = 0. Contributions (II) and (III) contain 
the coUinear divergence which starts at and (IV) is the overall infrared 
divergence with a pole. Therefore these have to be calculated in d = 4 — 2£ 
dimensions. 



4.2 Divergent parts 

Contribution (IV) is easy: it is simply a product of two one-loop diagrams, 
where the k loop is completely massless: 



iw = -^^'^'^T^(iriy^(-«' - ^o)^c(o,o,o) . (27) 

The expansion of the F functions starts with so C(0,0, 0) has to be cal- 

culated up to O(e^). The last order is done numerically. 

Contribution (II) is more complicated, since the I and k integrations do not 
decouple. However, the integrand is independent of z, the angle between l± 
and k±, so there is no function with a cut in the complex ki plane. We have 

d-2 +00 00 +00 

^11 = / dA:o(C(fco,0,0)-C(0,0,0)) J dt J dki^-^'^8) 

V 2 y —00 —00 

and perform the ki integration with Cauchy's theorem. We close the contour 
in the lower half plane, and there is only a contribution from P5 in the finite 
interval < ko < e, outside this interval all poles are on the same side of the 
real axis. 

The integration over t is of the form 

00 d-6 



so after a substitution ko = —ze the result is 

r(¥) 

1 

J dz^[z{l - z)]^ {C{ze,0,0) - C{0,0,0)) . (30) 


The expansion of the T functions starts at 1/e, so the integral over z, which 
will be done numerically, has to be evaluated up to 0{e). 
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Contribution (III) can be handled analogously by performing the /cq inte- 
gration with Cauchy's theorem and gives 





iri 




ri 


m 


1 



1 

dz^[z(l-z)]'^(C(0,ze,0)-C7(0,0,0)) . (31) 







For the symmetric case mi = m2, Im is equal to In. 



4.3 The finite part 

The most difficult is the finite part from contribution (I). In principle, the 
calculation is the same as in j^, |lO|, but care has to be taken due to the divergent 
behavior of the individual parts of this contribution. The integration over z, 
the angle between l± and k±, is elementary for each for the four terms. In fact, 
only the first term with the full C{ko, ki,t) dependence depends also on z, and 
therefore gives a result with a cut in the complex li and ki plane. 

The h and ki integrations are done with the aid of the residue theorem. 
Each term is still convergent on its own, so the contours can be closed indepen- 
dently. For the C{kQ, ki,t) term we have to take into account the cut which is in 
the upper half plane if /q + < and in the lower half plane if /q + > 0. By 
checking the half plane where the propagators have their poles, we find three 
contributing triangles in the {Iq, k^) plane from the pairs (Pi, P5), {P2, P^) and 
(P2,-P6), depicted in Fig. |. 

For the second term, C{kQ, 0, 0), we can close the contours arbitrarily. If we 
choose to close them in the same way as above (this is equivalent to avoiding 
poles from P3, the propagator that depends on both loop momenta) we find 
the same triangles from the same propagators as above. 

With C(0,A;i,0) we find something new in our parallel/orthogonal space 
technique: regardless how we close the contours, we end up with an unbounded 
area in the {lo,ko) plane. If we continue with the strategy of avoiding poles 
from P3, we have to close the h contour in the upper half plane if /q > and in 
the lower half plane if Iq < 0, furthermore we close the ki contour always in the 
upper half plane. This gives us the three areas shown in Fig. || from {P2,Pi), 
{P2,P5) and {P2,P6) respectively. 

For last term C(0, 0,0), which factorizes in / and k we can close so that in 
the h integration only P2 contributes, and in the ki integration only P5. This 
results in a square (Fig. |6|). 

Next we will do the t and later the s integration analytically. Since we left 
niQ finite, after partial fractioning the integrals are of the form 

00 00 

f f 1 1 1 

ds dt — tt: ,, = for term 1 

J J s + So - ^0 t + to - «0 ^ (at + b + i{) + csY - 4.st 
00 
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oo oo 

f 1 I 1 

ds / dt ; for term 2 & 3 

J s + So — iO t + to — iO at + b + iO + cs 



oo oo 

ds / dt ; for term 4. (32) 

J s + so-iOs + s'r,-iOt + to-iOt + t'r.-iO 



Several of these integrals become divergent as mg 0, whenever the cor- 
responding to is proportional to niQ. Consequently, the leading term is pro- 
portional to log(m6). However, as can be verified, pairs of the above integrals 
remain finite in the limit rriQ 0. The pairs depend on the position in the 
(Zq, ^o) plane (area A, B or C, Fig. ^. 



5 Numerical Results 

Tables 1-3 summarize the results. Generally, there are two sets of results 
for the momentum expansion method, namely two [n/n] Fade approximants 
involving 2n + l Taylor coefficients, which document that the results are stable 
and converge. It is the second column which should be compared with the 
numerical data. 

Let us now discuss the three cases separately. Case 7 with equal masses. 
Table 1, was the most challenging for the numerical method. Below the thresh- 
old, we could achieve sufficient numerical accuracy. Above the threshold, the 
numerical method could only achieve an accuracy of ~ 1%. Fortunately, for 
the equal mass case, M. Spira could kindly provide data which were used for 
comparison with the momentum expansion method to high accuracy. The lack 
of accuracy in the numerical method in this case is due to the degenerate cut 
structure which is present in the equal mass case. As a result, we are confronted 
with two large contributions which almost cancel. 

In Case 7 with different masses. Table 2, we chose the masses mi = 7712 = 
80GeV{^ My/) and mg = 180G'ey(~ Mt) and /i = mi. If we wish to calculate 
the process Z ^ bb the kinematics of interest is q'^/ml ~ 1.3, i.e. far below the 
threshold at /m\ = 4. In such a situation the momentum expansion method 
always seems to be superior to any other approach (to achieve 10 decimals 
precision with a few coefficients only is no problem). Nevertheless, also for low 

the numerical method yields high precision as well which will be sufficient 
for all practical purposes. Indeed the agreement is quite impressive. 

For higher q^ the momentum expansion method naturally gets less precise, 
in particular near the thresholds: here q^ = Ami ^^'^ 'f' — ("^1 + "^e)^) i-e- 
q^ /m\ ~ 10.6. Between the thresholds and above, both methods again show 
surprisingly good agreement. 

The situation is indeed quite different for Case 8. While for the case of 
interest in Z — > 66, i.e. q^ /m\ = 1 both approaches yield approximately the 
same precision, the momentum expansion method looses quickly for the higher 
q^, i.e. many more Taylor coefficients would be needed. 
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Quite generally speaking the results of the momentum expansion method 
obtained so far seem to indicate that the more heavy masses in a diagram 
occur, the better the convergence of the Pade approximants. This is at least 
qualitatively plausible since we start from a large mass expansion. It is also 
numerically verified comparing Case 7 and 8. 

In any case, splitting off systematically zero-mass thresholds yields more 
precise results for the momentum expansion method than keeping small but 
non-vanishing masses resulting in low thresholds. This has been discussed in 

for Case 5 (i.e. = rriQ = M and all others 0, see Tables 3 and 4 of 
that Ref.). For the numerical method to handle such cases is much easier. In 



fact, one could rely on the original integral representations 10 1 without the 
modifications considered here and achieve much better numerical results for all 
cases. 

Concerning the CPU-time, comparison of the methods is difficult. In the 
momentum expansion method the real job is to calculate the Taylor coefficients. 
Once they are known, for any (in the complex plane) the calculation of the 
diagram under consideration is a matter of seconds. Thus, providing Taylor co- 
efficients for the diagrams is something which practically settles the calculation 
of the diagrams under consideration once for all in a wide range of q"^. This 
is in particular true if only one nonzero mass is involved since in this case the 
Taylor coefficients are just rational numbers plus some well known irrationals 
like C(2), C(3) etc. The methods to calculate Taylor coefficients are at present 
still a matter of intense investigation (see e.g. also ||l^) and the hope is that 
we will soon have more efficient methods for their calculation and higher coef- 
ficients than used here can be obtained in the future. 



For the numerical method, we can give a rough estimate by comparison 
with the results of Q. Due to the extra subtraction terms, the CPU time 
needed here is about an order of magnitude larger than in the IR- and collinear 
convergent cases studied in 



6 Conclusion 

Two different methods to calculate scalar two-loop vertex diagrams with in- 
frared and collinear divergences have been investigated and interesting tech- 
niques for their calculation have been developed. An important case with dif- 
ferent masses has been solved, which is a step forward to a realistic calculation 
of the two-loop decay Z — > bb. The results show that both methods deliver con- 
sistent results. This confirms the expectation that the two techniques developed 
are applicable in the demanding cases considered here. 
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Table 1. The finite part for Case 7 with equal masses. 





fio/iol 




[14/141 

L / J 




numerical results 






Rc 


Im 


Rc 


Im 


Rc 


Im 


1.0 


2.764761286 


-0.1635351018 


2.764761286 


-0.1635351018 


2.76477 


-0.16353 


2.0 


1.619771096 


0.4534319337 


1.619771096 


0.4534319337 


1.61977 


0.45344 


3.0 


1.477131552 


0.5068403285 


1.477131552 


0.5068403285 


1.47713 


0.50684 


3.9 


3.0354569 


0.28920132 


3.035456900 


0.2892013177 


3.03546 


0.28920 


3.99 


5.1725 


0.0143 


5.17259 


0.014233 


5.17253 


0.014236 


4.01 


6.188 


2.748 


6.187 


2.749 


6.18778 


2.74948 


4.1 


3.27648 


3.92262 


3.27645 


3.92263 


3.27646 


3.92263 


5.0 


-0.438728 


2.54547 


-0.438728258 


2.545472409 


-0.438728 


2.54548 


6.0 


-0.89277747 


1.49214654 


-0.8927774793 


1.492146568 


-0.892777 


1.49215 


7.0 


-0.90571711 


0.92720655 


-0.9057171008 


0.9272065484 


-0.905717 


0.927207 


8.0 


-0.828158487 


0.59820497 


-0.8281584836 


0.5982049870 


-0.828159 


0.598205 


10.0 


-0.65092293 


0.26109287 


-0.6509229304 


0.2610928562 


-0.650923 


0.261093 


40.0 


-0.07133439 


-0.0467805 


-0.0713343755 


-0.0467805683 


-0.0713345 


-0.0467805 


100.0 


-0.0126312 


-0.0182863 


-0.012631500 


-0.018286325 


N/A 


N/A 


400.0 


-0.000556 


-0.002743 


-0.00055586 


-0.00274078 


-0.000555895 


-0.00274078 



Table 2. The finite part for Case 7 with diff'erent masses. 





[12/12] 




[15/15] 




numerical results 






Re 


Im 


Re 


Im 


Re 


Im 


1.0 


1.255976034 


0.4303855857 


1.255976034 


0.4303855857 


1.255974 


0.430385 


2.0 


0.6127661221 


0.4991658334 


0.6127661221 


0.4991658334 


0.61277 


0.499164 


3.0 


0.4346988621 


0.4502907526 


0.4346988621 


0.4502907526 


0.43470 


0.450292 


3.9 


0.4583945 


0.34378450 


0.458394486 


0.343784488 


0.4584 


0.34379 


3.99 


0.4926 


0.2501 


0.49254 


0.25004761 


0.4925 


0.25003 


4.01 


0.637 


0.190 


0.6359 


0.1902 


0.63598 


0.18994 


4.1 


0.73506 


0.43258 


0.73504 


0.43260 


0.73526 


0.43271 


5.0 


0.32211 


0.855621 


0.322110 


0.855622 


0.32219 


0.85569 


6.0 


0.00752 


0.802352 


0.0075200 


0.8023528 


0.00756 


0.80237 


7.0 


-0.16178 


0.682114 


-0.161777 


0.682110 


-0.16176 


0.68214 


8.0 


-0.25256 


0.56468 


-0.252547 


0.564693 


-0.25255 


0.56469 


10.0 


-0.3259 


0.3799 


-0.3256 


0.37987 


-0.32559 


0.37975 


40.0 


-0.06956 


-0.02984 


-0.0695667 


-0.029828 


-0.06956 


-0.02983 


100.0 


-0.01350 


-0.01529 


-0.013493 


-0.015279 


-0.01349 


-0.01528 


400.0 


-0.00072 


-0.00253 


-0.00071 


-0.002537 


-0.00071 


-0.00254 



Table 3. The finite part for Case 8. 



q^/m^ 


[12/12] 




[15/15] 




numerical results 






Re 


Im 


Re 


Im 


Re 


Im 


0.5 


81.17501719 


12.06458720 


81.17501719 


12.06458720 


81.1750 


12.0644 


1.0 


17.766 


19.97799834 


17.7659 


19.97799834 


17.7658 


19.9779 


2.0 


-0.605 


7.376 


-0.6047 


7.3759 


-0.604 


7.376 


3.0 


-1.753 


3.182 


-1.7543 


3.1815 


-1.754 


3.182 


4.0 


-1.596 


1.606 


-1.595 


1.603 


-1.595 


1.604 


5.0 


-1.33 


0.881 


-1.327 


0.880 


-1.326 


0.880 


6.0 


-1.10 


0.50 


-1.096 


0.503 


-1.095 


0.503 


7.0 


-0.92 


0.28 


-0.914 


0.289 


-0.913 


0.288 


8.0 


-0.77 


0.15 


-0.773 


0.159 


-0.771 


0.159 


10.0 


- 0.56 


0.01 


- 0.572 


0.023 


- 0.570 


0.025 
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Figure 2: Propagator numbering, external and loop momenta 
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Figure 5: Contributions from term 3 
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Figure 6: Contribution from term 4 
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Figure 7: Integration regions in the (/q, ko) plane 
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